Additional Boundary Conditions for Nonconnected Wire Media 
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Following our recent work [New J. Phys. 10, 053011, (2008)], here we demonstrate that due to 
strong nonlocal effects additional boundary conditions are essential to characterize the reflection 
of electromagnetic waves by nonconnected wire arrays using homogenization methods. Based on 
simple physical considerations, we derive the additional boundary conditions for the case where the 
wire medium is adjacent either to a dielectric or to a conducting material, and demonstrate that in 
the lossless case such boundary conditions ensure the conservation of the power flow. It is shown 
that the number of additional boundary conditions is related to the number of metallic wires in a 
unit cell. We illustrate the application of the novel boundary conditions to several configurations 
with practical interest. 

PACS numbers: 42.70.Qs, 78.20.Ci, 41.20.Jb, 78.66.Sq 



I. INTRODUCTION 

Artificial materials formed by long metallic wires have 
attracted significant attention in recent years [1-14]. In 
some circumstances such wire media may be character- 
ized by a plasmonic-type electric response, which poten- 
tially may enable interesting phenomena to occur such 
as superlensing [15], negative refraction [12, 13], or the 
excitation of strongly localized electromagnetic modes. 
However, it has been known for some time that in general 
wire media may fail to completely mimic the properties 
of a local plasma, and may have a spatially dispersive re- 
sponse [4, 6, 16-22], i.e. the electric displacement vector 
in a given point of space cannot be written exclusively 
in terms of the macroscopic field at the considered point, 
but, on the contrary, may depend on the distribution 
of the electric field in a neighborhood that encompasses 
many unit cells. Such property is a consequence of the 
fact that the metallic wires are spanned over many unit 
cells, and consequently the polarization acquired by the 
metallic inclusions in a given unit cell, (which is roughly 
proportional to the current along the wires), may de- 
pend significantly on the electric field outside the con- 
sidered cell. Recently, it was suggested that the spatial 
dispersion effects may be tamed either by attaching con- 
ducting structures on the wires or by coating the wires 
with a magnetic material [23]. In Ref. [15] it was also 
demonstrated that at infrared and optical frequencies the 
nonlocal effects may be significantly weakened when the 
radius of the metallic rods is comparable to (or smaller 
than) the skin depth of the metal. 

Even though in some applications spatial dispersion 
effects may be regarded as undesirable, in other circum- 
stances they may open new possibilities. For example, 
in part due to nonlocal effects, an array of parallel wires 
may behave as a material with extreme optical anisotropy 
and this may enable the transport and manipulation of 



'Electronic address: mario.silveirinha@co.it.pt 



the near- field [7, 9-11, 24, 25]. Moreover, in a recent se- 
ries of works it was shown that the spatially dispersive 
properties of a crossed wire mesh may enable the real- 
ization of materials with an extreme index of refraction 
[14], the realization of ultra-subwavelength waveguides 
[14, 26], superlensing [27], broadband all-angle negative 
refraction [28], and low-loss broadband anomalous dis- 
persion [29]. Some of these effects (e.g. the low-loss 
anomalous spectral dispersion) are specific of spatially 
dispersive materials, and cannot be observed in local ma- 
terials. Due to these and other potential applications, 
the accurate characterization and modeling of nonlocal 
materials gains increasing importance. 

One of the peculiarities of nonlocal materials is that 
they may support additional waves, i.e. for a fixed direc- 
tion of propagation and a fixed frequency, the number of 
plane waves supported by the material may be greater 
than two, differently from conventional local materials. 
For example, some of the plane waves may have identi- 
cal polarization (i.e. the same orientation for the electric 
field) and be associated with different wave vectors. It is 
well known that such property implies that the classical 
boundary conditions, which impose the continuity of the 
tangential components of the electric and magnetic fields, 
may be insufficient to solve a scattering problem at an 
interface between a spatially dispersive material and an- 
other material, due to the degrees of freedom associated 
with the "new" waves. A possible strategy to overcome 
this problem is the introduction of additional boundary 
conditions (ABC) [30-33]. In Refs. [5, 34] alternative 
"ABC-free" strategies are described. 

There is no general theory for the derivation of ABCs. 
In fact, the ABCs depend on the very specific microstruc- 
ture of the material, and should describe the dynamics of 
the internal variables of the nonlocal material, so that the 
degrees of freedom are removed. In recent works the ABC 
approach was successfully used to model the reflection 
and refraction of electromagnetic waves at an interface 
between a wire medium formed by parallel wires and a 
conventional (dielectric or conducting) material [35, 36]. 

The objective of this paper is to generalize the theory 
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of our previous works [35, 36] to the case of nonconnected 
"double" and "triple" wire media, and demonstrate that 
the introduction of suitable ABCs enables the accurate 
modeling of these structures using homogenization tech- 
niques. As mentioned before, such materials may have 
interesting potentials in several problems [14, 26-29], 
which justifies, besides obvious theoretical motivations, 
the present study. It should be mentioned that it is not 
trivial to extend the theory of Refs. [35, 36] to the case of 
more complex metamaterial topologies. For example, it 
will be shown here that in general at an interface between 
a crossed wire mesh and a regular dielectric two distinct 
ABCs are required, whereas for an array of parallel wires 
one single ABC is sufficient to characterize the interface 
effects [35, 36]. 

This paper is organized as follows. In section II, we 
derive the ABCs. In section III it is demonstrated that 
in the lossless case the proposed ABCs ensure the con- 
servation of the power flow through a wire medium slab. 
A formula for the Poynting vector in the spatially dis- 
persive material is derived. Then, in sections IV- VI, we 
illustrate the application of the ABCs in several scenar- 
ios of interest, and demonstrate their accuracy using full 
wave numerical simulations. In particular, we discuss 
the potentials of a wire mesh in the realization of ultra- 
subwavelength waveguides and the emergence of negative 
refraction. Finally, in section VII the conclusions are 
drawn. It is assumed that the fields arc monochromatic 
and have the time dependence e^*'^*. 

In this work, we use several different notations for the 
field entities. For the convenience of the reader, the 
meaning of these notations is summarized below for the 
electric field case: 



• e - microscopic electric field. This is the exact solu- 
tion of the Maxwell's equations, taking into account 
the exact microstructure of the material. 



• E - "bulk" electric field. This is the electric field 
that results from the homogenization of the bulk 
material. Specifically, the macroscopic field is ob- 
tained by averaging the microscopic field over a vol- 
umetric region within which the structure is invari- 
ant to translations along three independent direc- 
tions of space. 



• Eav.T - transverse averaged electric field. This defi- 
nition requires that the considered structure is peri- 
odic along two directions of space (assumed parallel 
to the xoy plane). However, the periodicity along 
z is not required. The transverse averaged field is 
obtained by averaging the microscopic field exclu- 
sively along the directions parallel to the xoy plane, 
as detailed in section II. 



II. DERIVATION OF THE ADDITIONAL 
BOUNDARY CONDITIONS 

Here, we generalize the ideas of [35, 36] and derive ad- 
ditional boundary conditions that enable the characteri- 
zation of the internal variables responsible for the spatial 
dispersion effects in nonconnected wire media. 




FIG. 1: (Color online) Geometry of the unit cell of the "bulk" 
double wire medium (A'' — 2). The cubic unit cell contains 
two nonconnected metallic wires. The dashed line represents 
an hypothetical cut of the bulk material with a plane normal 
to the z-direction. 

The material is formed by N mutually orthogonal and 
nonconnected sets of wires, oriented along the directions 
of space u„ {n = 1..N) [18], where N = 1,2,3 is the 
number of components of the wire medium. Each wire 
array is arranged in a square lattice with lattice constant 
a. The radius of the wires is r^,. For N = 1 the material 
consists of a single array of parallel wires, and for TV = 2 
{N = 3) it is formed by a double (triple) wire array. The 
distance between adjacent orthogonal wires is a/2. It 
is assumed that the wires are good conductors, i.e. the 
radius of the wires is much larger than the skin depth 
of metal at the frequency of operation. For reasons that 
will be clear shortly, in general the unit vectors u„ are 
not oriented along the coordinate axes x,y,z. 

As an example, the unit cell of the bulk double wire 
medium (crossed wire mesh) is represented in Fig. 1. 
The unbounded material is formed by the periodic repe- 
tition of the primitive cell. Since the wires intersect the 
boundaries of the unit cell, it should be clear that such 
construction effectively yields infinitely long wires. 

Let us suppose that the bulk material is sliced to form 
a planar interface with another adjacent material (a reg- 
ular dielectric or metal; see Fig. 4c). This operation 
breaks the translational symmetry of the system along 
the direction normal to the interface, which without loss 
of generality is assumed to be the z-axis. In particular, 
the wires lying in cells of the wire medium adjacent to 
the interface are sliced as well (except, possibly, if some 
of the wires are parallel to the interface). However, the 
periodicity of the system in the transverse plane {xoy- 
plane) is preserved. 
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Wc use the transverse-average (TA) field approach in- 
troduced in [37], in order to describe the electromagnetic 
wave propagation using homogenization methods. The 
TA-fields are obtained by averaging the "microscopic" 
fields over the directions parallel to the interface (x and 
y directions), i.e. over the directions along which the 
system has translational symmetry. It is assumed that 
the microscopic fields have the Floquet-Bloch property 
in the xoy plane, being characterized by the transverse 
wave vector k|| = {kx,ky,0). This situation occurs, for 
example, when the wire medium slab is illuminated by a 
plane wave, being in such circumstances k|| determined 
by the direction of incidence and by the frequency of op- 
eration. 

Following Ref. [37], the TA-electric and induction 
fields are defined as: 



Eav,T {z) 



Acell J 



e(r) e-"'^i ''dxdy 



Bav,T {z) = ^ [h (r) e-"'! i ''dxdy (1) 

^cell J 



where e and b are the microscopic electric and induc- 
tion fields, r^T represents the transverse unit cell in each 
z = const, plane, and A^eii is the area of I^t- As demon- 
strated in Ref. [37] , the TA-fields verify the the differen- 
tial system (assuming the time convention e~"^*): 



dz 
d 



Mo 



ik|| -I- — I X Eav,T = »wBav,T 



(2) 



where Sh is the relative permittivity of the host medium, 
and J(i,av is the averaged microscopic current. It should 
be clear that the TA-fields depend exclusively on z the 
coordinate. 

For the considered wire media, formed exclusively by 
metallic inclusions, the averaged microscopic current is 
given by [37]: 



1 



Jc (r) e 



-ikii 



1 



\U X Uz 



■dl (3) 



aA{z) 



where Jc = x b//io is the density of electric current on 
the surface of a given wire, v is the unit vector normal to 
the surface of the wire, OA [z) is the contour determined 
by the intersection of the surface of the wires in the unit 
cell and the pertinent z = const, plane, and dl is the el- 
ement of arc. Strictly speaking formula (3) is only valid 
when there are no wires parallel to the transverse plane 
(otherwise the integrand becomes singular), which is the 
case of interest in this work. The key result is that, simi- 
lar to the analysis of Ref. [36], it is possible to write Jd,av 
in terms of the microscopic electric currents induced on 
the metallic wires. 



Indeed, assuming that the wires are relatively thin, 
^ a, it is a good approximation to consider that 
the density of current over the n-th wire in the unit cell 
(parallel to the unit vector u„) is of the form (thin wire 
approximation): 



In{z) 



(4) 



where dDn represents the surface of the considered wire, 
/„ is the electric current, and n = 1..N, where A'^ is the 
number of components of the wire medium. Notice that 
the density of current is modulated by the exponential 
factor e+*'^ii '", because of the assumed Floquet period- 
icity of the microscopic fields in the transverse plane. 
Substituting Eq. (4) into Eq. (3), it may be shown that 
the averaged microscopic current is given by: 



(5) 



Thus Jd.av has indeed a simple relation with the micro- 
scopic electric currents. On the other hand, from Eq. (2) 
it is straightforward to write the averaged current as a 
function of the macroscopic TA-electric field: 

iw/ioJd.av [z) = 



ik\\ + — u 
" dz 



«k|| + —Viz I -Eav.T 



(6) 



Eqs. (5) and (6) establish the connection between the mi- 
croscopic currents and the macroscopic electric field and 
will be used in what follows to derive additional bound- 
ary conditions for nonconnected wire media. 

Let us consider first the case where the material adja- 
cent to the wire medium is non-conducting (e.g. a dielec- 
tric). In such situation, similar to the case of an array 
of parallel wires [35, 36], the microscopic electric current 
must vanish at the interface. This means that 7„ = 0, 
n = l..A^, at the interface with the dielectric material, 
i.e. the currents associated with different wires in the 
unit cell must vanish independently. But, Eq. (5) shows 
that this property implies that the averaged current ver- 
ifies: 



Jd,avU„ = 0, (diel. interface), n=l..N 



(7) 



This equation, along with Eq. (6), may be regarded as 
a set of A'' independent additional boundary conditions 
at the interface with the dielectric material. It should be 
clear that Eq. (7) may be written exclusively in terms of 
the macroscopic electric field and its derivatives, and that 
the number of additional boundary conditions is equal 
to the number of metallic wires in a unit cell {N = 2 
for the double wire medium and A/' = 3 for the triple 
wire medium). Obviously, Eq. (7) assumes that all the 
metallic wires intersect the interface. When some Un 
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is parallel to the interface, in general J^.av-Un does not 
vanish at the interface because the path of the current is 
not interrupted. 

The ABCs (7) are a generalization of the results of 
our previous works [35, 36] for an array of parallel wires 
(A^ = !)• As demonstrated in [35, 36], when = 1 the 
ABC and the classical boundary conditions imply the 
continuity of the normal component of the electric field 
multiplied by the host permittivity at the interface. In 
particular, when the host material is air, the ABC, to- 
gether with the classical boundary conditions, is equiva- 
lent to the continuity of all the Cartesian components of 
the electromagnetic field. Interestingly, when N = 2 oi 
N = 3 the situation is quite diff'erent, specifically, since 
the number of ABCs is greater than one, the continuity 
of the normal component of the electric field multiplied 
by the host permittivity (one single independent equa- 
tion) cannot be regarded as equivalent to all the ABCs 
(two or three independent equations). 

Another case of interest occurs when the material adja- 
cent to the wire medium has very high conductivity (e.g. 
a metal). It is obvious that in such configuration the 
ABCs (7) arc not valid, because, assiiming that the wires 
are connected with good Ohmic contact to the ground 
plane, the current path is not interrupted at the inter- 
face. In fact, it was demonstrated in Ref. [36] that it is 
the microscopic electric density of charge Cc that vanishes 
at the connection point with the metallic surface, rather 
than the electric current. Thus, electric charge cannot be 
accumulated at the connection points between the wires 
and the ground plane. For a wire directed along the n-th 
direction, the condition Cc = is equivalent to 



where Jc 



inMp+ikii.r 



= 0, 



is the density of current along 



the considered wire (see Eq. (4)), and s is a coordinate 
measured along the wire axis. It is clear that for the n- 
th wire ^ = u„.V, and hence the condition = 

is equivalent to u„. (ik|| + u^^) /„ (z) = 0. Thus, us- 
ing Eq. (5), it is readily found that the averaged current 
must verify: 



ik|i + Uj 



d 
dz 



.U„U„.J 



n »-*n.«' ti.av 



0, (metallic interface) 
n=l..N (8) 



The above equations define a set of N independent ad- 
ditional boundary conditions at the interface between a 
nonconnected wire medium and a ground plane, and are 
a generalization of the results of [36] . Again, it should be 
clear that Eqs. (6) and (8) define a functional relation 
between the macroscopic electric field and its derivatives 
at the interface, and that it is implicit that all the wires 
intersect the interface. 

The application of the proposed ABCs will be illus- 
trated in sections V and VI, where, for the case of a 
double wire medium, it is shown with full wave simula- 
tions that the ABCs may enable the accurate analytical 
modeling of wave propagation. 



An important point, which is discussed next, is the 
characterization the TA- macroscopic fields in the wire 
medium. The obvious idea is to write the macroscopic 
fields in the material as a superposition of plane waves. 
The plane waves may be characterized using the homog- 
enization model proposed in [14, 18, 19]. Specifically, the 
dielectric function of the nonconnected wire medium is, 

N 

f(w,k) = e/jl-h^(e„,„(w,k)-e/j)u„u„ 



en,n (W,k) = Eh 1-1- 



fv{Sm/Sh-i) 



(9) 



where I is the identity dyadic, u„u„ = u„ u„ 
represents the dyadic (tensor) product of two vectors, 

(jp = [27r/(ln(a/27rr^) -F0.5275)]^/V« is the plasma 
wavenumber, fy = -k [r^/af, Sh is the host relative 
permittivity, Sm = £m {^) is the metal relative com- 
plex permittivity, c is the speed of light in vacuum, 
k = (fcx) ky, kz) is the wave vector, and fc„ = k.u„. 

Using the dielectric fimction (9) it is possible to cal- 
culate the bulk macroscopic fields (E, B). An impor- 
tant observation is that, for a general metamaterial, the 
bulk macroscopic fields may not be coincident with the 
TA-fields used in the formulation of the ABCs (7) and 
(8). Indeed, the bulk macroscopic fields, as defined in 
[18, 37], are obtained by averaging the microscopic fields 
over the unit cell (a volumetric region) of the periodic 
material, whereas the TA-ficlds are obtained by averag- 
ing the microscopic fields over the transverse unit cell (a 
surface) [37]. Fortunately, as demonstrated in Appendix 
A, for nonconnected wire media, and supposing that all 
the wires intersect the interface (i.e. there are no wires 
parallel to the interface), it is possible to identify to a very 
good approximation the bulk macroscopic fields with the 
TA-fields. This means that it is possible to drop the sub- 
scripts "av,T" in Eq. (6), and characterize the TA-field 
Eav,T using the dielectric function of the bulk material. 
These ideas are further clarified in sections V and VI. 



III. CONSERVATION OF THE POWER FLOW 

An important question that may be formulated is if in 
the lossless case (i.e. when both Sh and Sm are real val- 
ued) the proposed ABCs ensure the conservation of the 
power fiow in the wire medium. Next, we demonstrate 
that that is indeed the case. The theory developed below 
is based on the bulk medium fields, which as mentioned 
in the end of the previous section, may be regarded as 
equivalent to the transverse averaged fields in the case of 
interest. 

Let (E, H) be a solution of Maxwell's equations in the 

homogenized bulk wire medium such that the variation 
in the x and y coordinates is of the form e**^^! "", with 
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k|| = (kx,ky,0) a real vector. The bulk magnetic field is 
by definition H = B/yUg- Thus, (E,H) verify: 



V X H = -iuD 

V X E = iw/ioH 



(10) 



where D = £oe/iE + P is the electric displacement vector 
in the spatially dispersive material and P is the nonlo- 
cal (generalized) polarization vector (relative to the host 
material). Let us define S° as follows: 



1 



Re{E X H*} 



(11) 



It is well-known that in general S" cannot be identified 
with the Poynting vector in a spatially dispersive mate- 
rial [38, 39]. In fact, there is no general formula for the 
Poynting vector in a spatially dispersive material. The 
only case for which the Poynting vector is actually known 
is when the material is lossless, and the electromagnetic 
field is associated with a propagating plane wave [39] . In 
what follows, we will demonstrate that for lossless non- 
connected wire media, it is possible to define unambigu- 
ously the 2;-component of the Poynting vector, even if 
the electromagnetic field has an arbitrary variation with 
z (not necessarily a propagating plane- wave). 

To begin with, we note that from Eq. (10) it is imme- 
diate that: 



1 



1, 



V.S" = -Re {-zwE.D*} = -Re {-iwE.P*} (12) 

where the second identity is a consequence of the fact 
that the host permittivity en is a real number in the 
absence of loss. The next step is to relate P and E. In 
the spectral (Fourier) domain we clearly have that: 



'(k)=eo(e(a;,k)-£/,l).E(k) 



(13) 



where e is the relative permittivity of the wire medium, 
the symbol denotes the spatial Fourier transform, 
and here k represents the Fourier coordinates. In par- 
ticular, the projections of E and P onto the princi- 
pal direction u„ of the dielectric function verify P„ = 
£:o{£n,n{<^,kn) - Sh) En with kn = k.u„ (n = 1..N). 
Thus, using Eq. (9) it follows that: 



En 



1 



£o£h(3p 



01 



fv {£m/£h - 1) 



Pn (14) 



Hence, calculating the inverse Fourier transform, we ob- 
tain the following relation in the space domain: 



En 



1 



£o£hPn 



-ZU„.V) - -^£h + J 7T 

fv {£m/£h - 1) 



Pn 

(15) 

But since, it is assumed that the variation of the electro- 
magnetic field in the x and y coordinates is of the form 



g'''!! '' we have that — iu„.V = u„.k 



after some simplifications, it is found that: 
1 



iu„.U2^. Thus, 



Re{-iP*En} 



£o£hPp 



,,Re^z(u„.u,)^P:^ 



n 



-2(k||.u„) (u„.u^)P; 



dP„ 
dz 



(16) 



Next we note that Re {iP„*i0^} = Re (iP*^^^)} 
and that ^e{P*^} = £ to finally write: 



Re{-iP„*^„} = 



^-Re|(u„.u.) ^P„ — 



(k||.U„) (Un-Uz) 



(17) 



We are now in a position to calculate the divergence 
of S°, given by Eq. (12). On one hand, we note that 

due to the assumed dependence of (E, H) on the x and 
y coordinates, it is clear that S'' depends exclusively on 

z, and thus V.S'^ = where 5° is the z component 
of S°. On the other hand, it is obvious that P*.E = 

N 

PnEn because the only non- vanishing components of 

n=l 

P are precisely the P„'s (remember also that the unit 
vectors u„ are mutually orthogonal). Therefore Eqs. (12) 
and (17) show that. 



where Sz is given by: 



dz 
1 







(18) 



2 £o£hPv 



N 



Re 



n=l 



dz 



•UrrPn 



(19) 



As discussed below Sz is the z component of the Poynting 
vector in the spatially dispersive material. Before that 
discussion, let us show that the conservation law (18) 
implies that the ABCs introduced in this work guaran- 
tee the conservation of power flow at an interface with a 
dielectric or a metallic material. 

Indeed, it is clear from Eq. (19) that at the points z 
such that for every n = l..A^, 

either P„ = or (^ik,, + -^u^^ .u„P„ = (20) 

we have that Sz = 5'°. Comparing Eqs. (2) with Eqs. 
(10), and using the fact that in the nonconnected wire 
medium the bulk medium fields can be identified with 
the TA- fields (see Appendix A), it is evident that the 
polarization vector verifies. 



P = 



1 



-ILO 



d,ave " 



(21) 
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where Jd,av is the averaged microscopic current. But 
then it follows that the first condition in Eq. (20) is 
equivalent to the ABC (7), whereas the second condition 
is equivalent to the ABC (8). Hence, we conclude that 
at the points z of the material where either the ABCs 

(7) (associated with dielectric interfaces) or the ABCs 

(8) (associated with metallic interfaces) are enforced, we 
have that Sy_=Sl. 

Let's now consider a truncated wire medium (wire 
medium slab) with interfaces at z = Zi and z — zf. 
From Eq.(18), in the lossless case Sz is constant inside 
the wire medium, and consequently when the ABCs are 
enforced, 5° has the same value at the two interfaces: 
Sz{zi) = S^{zf) (calculated from the wire medium side). 
But due to the classical boundary conditions (continuity 
of the tangential components of E and H at the inter- 
faces), may be evaluated either at the wire medium 
side of the interface or at the exterior side, being the 
result the same. Thus, we have demonstrated that the 
proposed ABCs together with the classical boundary con- 
ditions imply that S^{zi) = S^izf), being evaluated at 
the exterior side of the interface, where it is obviously the 
z component of the Poynting vector. This shows that in 
the lossless case the ABCs (7) and (8) (or more generally 
the condition (20)) guarantee, in fact, the conservation 
of the power flow through a wire medium slab, as we 
wanted to prove. This result is valid for simple (TV = 1), 
double {N — 2) or triple {N — 3) wire media. 

To show that Sz can be identified with the z compo- 
nent of the Poynting vector in the spatially dispersive 
material, we consider the typical case where the electro- 
magnetic field is a superposition of plane waves: 

E(r) = ^E,e+* (22) 
I 

where E; is a constant vector (determines the polariza- 
tion of the plane wave), and k; = k|| -f kz^Uz is the wave 
vector associated with the plane wave. The z-component 
fci'^ may be either real (positive or negative) or complex 
valued (as discussed in section IV, the wire medium may 
support evanescent modes even in the absence of loss). 
It is supposed that each individual plane wave verifies 
Eqs. (10). In Appendix B, we demonstrate that for such 
superposition of plane waves, Sz, defined by (19), may 
be rewritten as: 

Sz= J2 Re|(iEzxH;,].u, 

where the summation is restricted to the indices I , m such 
that k; = kj„, and H; is the magnetic field associated 
with E; . But for a single plane wave with real wave vector 
the above formula reduces to the well-known expression 
of the Poynting vector in a general spatially dispersive 



material [38, 39]. This demonstrates that 5*2 may in- 
deed be identified with the z-component of the Poynting 
vector. It should however be stressed that the above for- 
mula is more general than the results reported in [38, 39], 
which apply exclusively to a single plane wave with real 
wave vector. The result (23) generalizes the classical re- 
sult to the case of a superposition of plane waves (pos- 
sibly associated with complex wave vectors) in a lossless 
nonconnected wire medium. 



IV. CROSSED WIRE MESH 

In the rest of the paper we apply the developed theory 
to the particular case of a crossed wire mesh formed by 
a double array of metallic wires [6, 14, 19]. It is assumed 
that the wires are parallel to the a::oz-plane, and that 
ui = (l,0,l)/\/2 and Ua = (-1,0,1)/V2. Thus, the 
angle of the wire axes with respect to the interface normal 
is ±45°. A cut of the considered structure in the planes 
yoz and xoz is represented in Fig. 2. We will analyze the 
scattering and guiding of electromagnetic waves by the 
wire medium slab in the cases where the wave vector is 
confined to one of these two planes. Moreover, in case of 
propagation in the yoz plane with will assume that the 
electric field is polarized along the x-direction, whereas 
in case of propagation in the xoz-plane we will restrict 
our attention to the situation where the magnetic field 
is along the y-direction (thus, in both cases, the electric 
field is in the xoz plane). 




-3-2-10123 -3-2-10123 



k^a k^a 

FIG. 2: (Color online) Isofrequency contours of the funda- 
mental plane wave mode for (a) propagation in the yoz plane 
with electric field along the a;-direction; (b) propagation in 
the xoz plane with electric field in the same plane. The ra- 
dius of the wires is r„ = 0.05a and the metal is assumed 
perfectly conducting. The text insets indicate the value of 
the normalized frequency uaa/c. 

Before addressing the problem of propagation in a fi- 
nite structure (wire medium slab), next we will briefly re- 
view some key properties of the electromagnetic modes in 
the unbounded crossed wire mesh. Unlike a conventional 
local plasma, the metamaterial may support propagating 
modes for arbitrarily low frequencies [6, 8, 14, 18, 19]. 
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Specifically, in the lossless case (e.g. for perfect elec- 
trical conductors (PEC), i.e. = — cxd), and for fre- 
quencies much lower than the effective plasma frequency, 
w/c <C /3p/v^) the wire medium supports a propagating 
plane wave with electric field in the xoz plane. 

Following Refs. [14, 27], the plane waves supported 
by the material in case of propagation in the yoz plane 
{kx = 0) with electric field along the a;-direction are char- 
acterized by the dispersion characteristic: 

e {u},kz) (co/c)'^ = ky + kl, with 



e (w, k^) = eh [l 



/v(em/£h-l) 



1^1 



(24) 



The isofrequency contours of the propagating mode are 
depicted in Fig. 2a, for PEC wires with normalized ra- 
dius Tw/a = 0.05. Clearly, the material is strongly 
anisotropic, and the isofrequency contours remind el- 
lipses with a large axial ratio. It is seen that the con- 
tours are nearly perpendicular to the z-direction, and 
consequently the energy tends to flow along z. This sug- 
gests that adjacent wire planes tend to guide the wave 
along z, and obstruct the propagation along y. It was 
shown in Ref. [14] that the effective index of refraction 
for propagation along z is (assuming £/j = 1 and PEC 
wires): 



kzC 



fe„=0 




— (25) 



Since /3p ^ 1/a, one interesting feature of this structure 

is that for a fixed frequency the index of refraction may 
be greatly enhanced by reducing the spacing between 
the wires (maintaining the metal volume fraction) [14]. 
Other mechanisms to design materials with a high-index 
were proposed in [40, 41], but the required microstruc- 
turcs arc arguably much more complicated to fabricate 
than a crossed wire mesh. 

As could be expected from the geometry of the mate- 
rial, the propagation properties in the xoz plane (ky = 0) 
are very different. The dispersion characteristic is now 
[28]: 



+ 



^2 



- (w/c)^ £ii fc2 - (w/c)^ £22 



= 1 



(26) 



where £„,„ (n = 1, 2) is given by Eq. (9), k"^ = fc^ 
and kn = k.u„, with k = {kx,0,kz). The corresponding 
isofrequency contours are depicted in Fig. 2b, and remind 
two perpendicular hyperbolas. The electric field associ- 
ated with each plane wave is nearly (but not exactly) 
tangent to the isofrequency contours [28]. The shape of 
the isofrequency contours suggests that propagation is 
favored along the coordinates axes, i.e. when the wave 
vector makes an angle of ±45° with the wire axes. On 



the other hand, when the wave vector is normal to one of 
wire axes the structure has a directional bandgap. The 
hyperbolic shaped isofrequency contours may enable the 
emergence of negative refraction [6, 19, 28], as will be 
further discussed in section VI. 



V. PROPAGATION IN THE YOZ PLANE 

In what follows, we apply the ABCs derived in section 
II to characterize the scattering of plane waves by a dou- 
ble wire medium slab, assuming propagation in the yoz 
plane (Fig. 2a). As in section II, it is assumed that the 
normal to the interfaces is along the z-direction. 

To begin with, we calculate the averaged microscopic 
current in the wire medium slab. It is assumed that the 
excitation is such that the electric field is along the x- 
direction. Clearly, for propagation in the yoz-plane the 
transverse wave vector is such that k|| = kyUy, where ky 
is determined by the excitation (for an incoming plane 
wave, ky = oj/csinOi where 9i is the angle of incidence). 
Thus, from Eq. (6), it is found that the averaged current 
verifies, 



iu)IIo3d,a.v {z) 



(27) 



where is the electric field in the wire medium slab. 

The ABCs at a dielectric interface can now be easily ob- 
tained using Eq. (7). As discussed in section II, in the 
double wire medium case, one needs to impose two dis- 
tinct ABCs at the interface (A^ = 2). However, in the 
present case the two ABCs (7) impose the same condi- 
tion on the macroscopic electric field: 



(PEx 

dz^ 



+ (^£h'^- klj Ex = (diel. interface) (28) 



This degeneracy occurs due to the exceptionally high- 
symmetry of the system. In the general case, the ABCs 
are not redundant and yield, in fact, independent equa- 
tions, as will be shown in Sect. VI. Similarly, at a metal- 
lic interface, the two ABCs (8) also yield a single equa- 
tion: 



dz3 + 1 - 



dEx 
dz 



(metallic interface) 



(29) 

The equations (28) and (29) correspond exactly to the 
additional boundary conditions that were enunciated in 
our previous works [14, 27] (with no proof). These ABCs 
can be readily applied to solve a scattering or wave guid- 
ing problem when the wave vector is confined to the yoz 
plane. Since the use of the ABCs for this specific con- 
figuration was already addressed in Refs. [14, 27], our 
discussion here will be somehow abbreviated. 
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A. Slab standing in free-space 

First we consider the case where the wire medium 
slab stands in free-space, and is iUuminated by a plane 
wave characterized by the transverse wave number ky = 
uj/csmOi (see the inset of Fig. 3). As mentioned in 
Sect. II, the key assumption is that the electromagnetic 
field inside the wire medium slab can be written as su- 
perposition of plane waves characterized by the dielectric 
function of the bulk (unbounded) material. Because of 
translational invariance ky must be conserved, and thus 
the pertinent plane waves can be found from the solution 
of the dispersion characteristic (24) with respect to kz, 
with bj and ky determined by the excitation. As discussed 
in Refs. [14, 27], the solution of Eq. (24) yields exactly 
two different solutions for k^. It should be clear that for 
a regular local dielectric slab there is only one single solu- 
tion. This clarifies why for propagation in the yoz plane, 
a single ABC is sufficient to characterize interface effects 
(despite that the wire medium has two components), i.e. 
since there is only one additional wave only one ABC is 
required. Thus, the electric field inside the the metama- 
terial slab may be written as (the y-dependence of the 
field is suppressed): 

^^-g^»fe.,(2)^ (30) 

where kz^(i) — k^^^i) (w, ky) are the solutions of the dis- 
persion equation (24). As discussed in Refs. [14, 27], for 
long wavelengths and in the absence of loss, there is only 
one propagating mode in the wire medium, i.e. one of 
the propagation constants, let's say kz.(i), is real valued 
and is associated with the isofrequency contours depicted 
in Fig. 2a, while the other propagation constant, fcz,(2)j 
is pure imaginary. 

The electric field in the air regions (below and above 
the metamaterial slab) can be written as well as a su- 
perposition of (free-space) plane waves. By enforcing the 
continuity of the tangential components of the electric 
and magnetic field, and the ABC (28) it is possible to 
determine the reflection and transmission coefficients for 
plane wave incidence. For further details the reader is 
referred to Ref. [27] [in particular, formula (5) of Ref. 
[27] gives a closed form expression for the transmission 
coefficient T]. 

In the first example, we consider a metamaterial slab 
formed by PEC wires, and such that the spacing between 
the wires is a = L/15, where L is the thickness of the slab. 
The host material is air. The amplitude of the transmis- 
sion coefficient for plane wave incidence along 6i = 0.1° 
is plotted in Fig. 3 as a function of the normalized fre- 
quency. It is seen that the results obtained with the 
nonlocal homogenization model and the proposed ABCs 
(solid line) compare well with the results obtained with 
the full wave electromagnetic simulator CST Microwave 
Studio (discrete symbols) [42]. 

As mentioned before, in the absence of loss, the prop- 
agation constant k^^^i) may be assumed real valued (as- 




Local model 
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FIG. 3: (Color online) Magnitude of the transmission coeffi- 
cient as a function of frequency for a metamaterial slab with 
a lattice constant a = L/15 and = 0.05a and a fixed 
thickness L. The solid line represents the spatially dispersive 
(nonlocal) model, the dashed line represents the local model, 
and the discrete symbols were calculated with CST Microwave 
Studio. The inset shows the geometry of the problem. 

sociated with a propagating mode in the wire medium), 
whereas the propagation constant of fcz,(2) is pure imag- 
inary (associated with an evanescent mode). Thus, it 
might be thought that this later evanescent mode would 
play a minor role in the response of the metamaterial 
slab. In order to test this hypothesis, we have calculated 
the transmission coefficient obtained by setting = 
in Eq. (30). This condition removes the extra degree 
of freedom of the problem, and thus to calculate the 
transmission coefficient it is sufficient to impose the clas- 
sical boundary conditions at the interfaces, i.e. if one 
neglects the effect of the evanescent wave it is not nec- 
essary to consider an ABC at the interfaces. We refer 
to the results obtained with this approximation, as the 
"local model" results. A bit surprisingly, as reported in 
Fig. 3, the results obtained with the local model (dashed 
line) are qualitatively very different from the results ob- 
tained with the nonlocal model. In particular, the fre- 
quencies for which the local model predicts the maxima 
of T correspond to the minima of the nonlocal model, 
and vice-versa. This demonstrates that the effect of the 
evanescent mode and the ABCs cannot be neglected, and 
are essential to describe accurately the response of the 
material. It should be noted that the evanescent mode, 
even though strongly attenuated in the interior of the 
metamaterial slab, may be excited in the close vicinity 
of the interfaces in order to ensure that the microscopic 
current in the metallic wires vanishes at the interfaces. 
Thus, even though the wave propagation in the interior of 
the slab is dominated by the propagating mode, the "ef- 
fective wave impedance" at the interfaces is determined 
by both the propagating and the evanescent modes. This 
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explains the disagreement between the local and nonlocal 
models reported in Fig. 3. 

One of the interesting properties of the crossed wire 
mesh is that it may interact very strongly with the incom- 
ing wave, even when the the length of the metallic wires is 
electrically short, effectively behaving as a material with 
very large dielectric constant [27]. For instance, in the 
example depicted in Fig. 3 the first dip of the transmis- 
sion coefficient occurs for ujL/c « 0.2, which corresponds 
to the metallic wires with length = ^/2L = 0.04Ao. 
This value is one order of magnitude smaller than the tra- 
ditional Ao/2 resonance of a metallic wire. Moreover, as 
discussed in Ref. [27] , the resonance length of the wires 
may be made arbitrarily small by increasing the density 
of wires (number of wires per unit of volume), keeping the 
metal volume fraction unchanged. Indeed, unlike conven- 
tional metamaterial designs, the effective index of refrac- 
tion of the material does not saturate when the inclusions 
are scaled and the lattice constant a is made smaller. 
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FIG. 4: (Color online) Panel (a): generic metamaterial slab 
formed by inclusions whose size is smaller than the lattice 
constant. Panel (c): crossed wire mesh slab formed by inclu- 
sions that are spanned over many unit cells. As discussed in 
the main text, unlike a standard metamaterial, the response 
of the crossed wire mesh does not saturate when the density 
of inclusions (number of inclusions per unit of volume) is in- 
creased, keeping the volume fraction constant (panels (b) and 
(d)). 



To further clarify these ideas we consider the scenarios 
depicted in Fig. 4. In panels (a) and (b) we represent 
a conventional metamaterial slab (e.g. an array of metal 
or dielectric spheres) with fixed thickness L, but different 
lattice constants ax- The volume fraction of the inclu- 
sions is assumed to be the same in the two cases. Each 
inclusion is entirely contained in a basic cell. The in- 
coming wave illuminates the slab along the normal direc- 
tion. From the classical Clausius-Mossotti formula [43] , it 
should be clear that in the quasi-static limit the value of 
the index of refraction of the bulk metamaterial will even- 
tually saturate and become independent of the lattice 
constant a^, when ax is made sufficiently small (keeping 
the volume fraction of the inclusions unchanged). In fact, 
for a fixed frequency, and for ax increasingly small, the 
index of refraction will obviously converge to the static 
case value. In particular, it is expected that in the quasi- 
static limit the transmission coefficient in configurations 



(a) and (b) becomes independent of a^;, provided ax is 
sufficiently small. 

The situation for a crossed wire mesh is completely 
different. As discussed in Sect. IV - see Eq. (25) - for 
a fixed frequency and a fixed metal volume fraction, the 
index of refraction of the bulk wire mesh can be made ar- 
bitrarily large by reducing the lattice constant ax- Thus, 
despite the fact that the volume fraction of the metal is 
kept constant, the response of the structures depicted in 
Fig. 4c and Fig. 4d may be radically different. In partic- 
ular, unlike in a conventional metamaterial, the response 
does not saturate when ax is made increasingly small. 
This important qualitative difference between the wire 
mesh and a conventional metamaterial stems from the 
fact that in a wire mesh each inclusion is spanned over 
many unit cells, i.e. each individual wire crosses several 
unit cells (see Figs. 4c and 4d). Obviously, such prop- 
erty precludes the application of the Clausius-Mossotti 
formula, and clarifies the strong coupling between the 
different inclusions and the anomalous electric response 
of the wire mesh. 

To illustrate how the response of the wire mesh is mod- 
ified when the density of wires is increased, we have cal- 
culated the dispersion characteristic of the guided modes 
supported by the metamaterial slab (i.e. the waves that 
may propagate along the y-direction attached to the 
interfaces, even in the absence of an external source). 
The dispersion characteristic of the guided modes, ky — 
ky (w), may be obtained from the poles of the transmis- 
sion coefficient for plane wave incidence [27]. 




FIG. 5: (Color online) Dispersion characteristic of surface 
waves for a metamaterial slab with thickness L and different 
values of the lattice constant a. The radius of the wires is 
r„ — 0.05a. Solid lines: metamaterial slab stands in free- 
space. Dashed lines: grounded metamaterial slab. 

In Fig. 5 we depict (solid lines) the calculated disper- 
sion characteristics for a slab with thickness L and dif- 
ferent values of the lattice constant a. The metal volume 
fraction is constant in all the examples (r^, = 0.05a). The 
results were obtained using the homogenization model. 
Similar to a conventional dielectric substrate, the crossed 
wire mesh supports a guided mode with no cut-off fre- 
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quency, i.e. for arbitrarily long wavelengths. Consistent 
with the previous discussions, it can be seen that for 
a fixed frequency the efiiective index of refraction of the 
guided mode, kyc/to, increases significantly when the lat- 
tice constant a is reduced from a = L/S to a = L/15, i.e. 
when the density of the wires is increased. This simple 
example illustrates how unlike a conventional metama- 
terial slab, the response of the crossed wire mesh does 
not saturate as the lattice constant is made increasingly 
smaller. 

Even though, for simplicity, the examples considered 
here assume PEC wires, the proposed homogenization 
model can be applied as well when either the effect of 
loss or the plasmonic response of the metal are taken 
into account. For further details, the reader is referred 
to Ref. [27] , where the model has been validated against 
full wave simulations in such scenarios. 



B. Grounded Slab 



Next, we consider the case in which the bottom face 
of the metamaterial slab is covered with a metallic plane 
[14]. It is supposed that the metallic wires are attached 
to the ground plane with good Ohmic contact. In this 
configuration, assuming plane wave incidence with elec- 
tric field polarized along the x-direction, the electric field 
inside the metamaterial slab can still be written as in Eq. 
(30). The boundary conditions at the air interface are the 
same as in Sect. V A. On the other hand, at the metallic 
interface Ex must vanish, and the additional boundary 
condition (29) must be enforced. Proceeding in this man- 
ner (see Ref. [14]), it may be proven that the reflection 
coefficient referred to the air interface is: 



270 ikl 
1 + ^ 



(2) 



A-2 



D {uj,ky) 



kz,{2) ill + ^^,(2)) tan {kz,{i)L) - {jI + klj^^^ tan (A;^,(2)Z/) 



(31) 



D (w, ky) ^ A:^,(i)A:^,(2) l-yl (7^ + fc^ + fc^ (3^^ + fcf + 

( 



ll + kln)) ill + kl{2) 



^2,(1) + ^^,(2)) (^^.(1)-^) (^^,(2)-^) - 2A:^^(i)A:^^(2) sec (fc^,(i)L) sec (A:«,(2)-^) 



(32) 



70 {k%(2) - kl^i)) (7ft + fc^,(2)) tan (fcz,(i)i) - A;^,(i) (7^ + tan (^^,(2)-^) 

In Eqs. (31)-(32), fc^ ^i) and A:z,(2) are defined as in Sect. VA, L is the thickness of the slab, 70 = ^ky — u'^eqUo, 
7/1 = ^Jky — Lu'^eoHo^h, and ky is the transverse wave number of the incoming wave. 



To illustrate the application of the formula, we have cal- 
culated the phase of the reflection coeflicient as a function 
of frequency for a grounded metamaterial slab such that 
a = L/10 and different incidence angles. The wires are 
assumed perfectly conducting and thus the magnitude of 
p is unity. The results obtained with Eq. (31) are de- 
picted in Fig. 6 (solid lines) superposed on data obtained 
using CST Microwave Studio (discrete symbols) [42] . An 
excellent agreement is revealed supporting the validity of 
the proposed homogenization methods. It is interesting 
to note that the frequencies where the phase of p vanishes 
are nearly independent of the angle of incidence. Thus, 
at such frequencies the grounded slab may mimic very 
closely the behavior of an ideal perfect magnetic conduc- 
tor [44]. Moreover, due to the large index of refraction of 
the wire mesh such behavior may be obtained with a very 
subwavelength slab. For example, the first resonance in 
the example of Fig. 6 occurs for L = 0.02Ao. 

We have also calculated the dispersion characteris- 
tic ky = ky (w) of the guided modes supported by the 

groimded slab. The dispersion characteristic is obtained 
by numerically solving the equation D{ijj,ky) = with 



respect to ky, where D{u),ky) is defined by Eq. (32). 
The calculated results for a slab with fixed thickness L 
and different values of the lattice constant are depicted 
in Fig. 5 (dashed lines). Consistent with the results re- 
ported in Ref. [14] and the discussion of Sect. VA, it 
is seen that when the density of wires is increased (i.e. 
a is reduced) the guided mode becomes more attached 
to the metamaterial slab. Unlike the case where the slab 
stands in free-space the guided modes can only propagate 
above a certain cut-off frequency. Despite that, due to 
the anomalously high index of refraction of the crossed 
wire mesh, the metamaterial slab thickness may be very 
subwavelength at the onset of propagation of the funda- 
mental mode [14, 26]. For example, for a = L/10 the 
thickness of the slab at the cut-off frequency of the fun- 
damental mode is as small as L = 0.02Ao. Thus, the 
proposed structure may enable the realization of ultra- 
subwavelength waveguides, as demonstrated experimen- 
tally in [26]. 

It should be mentioned that, unlike a conventional di- 
electric slab, the dispersion characteristic of a grounded 
crossed wires slab cannot be obtained from the dispersion 
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FIG. 6: (Color online) Phase of the reflection coefficient as 
a function of the normalized frequency for different angles of 
incidence. The lattice constant is a = L/10 and the radius 
of the wires is — 0.05a, where L is the thickness of the 
grounded slab. The host material is air. Solid lines: Homog- 
enization model. Star and triangle shaped symbols: full wave 
simulations for d = 15° and 9 = 85°, respectively. 
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FIG. 7: (Color online) Amplitude of the reflection coefficient 
as a function of the transverse component of the wave vector 
ky, and different frequencies of operation ui. The lattice con- 
stant is a = L/10 and the radius of the wires is = 0.05a, 
where L is the thickness of the grounded slab. Solid lines: 
Homogenization model. Discrete symbols: full wave simula- 
tions. 



VI. PROPAGATION IN THE XOZ PLANE 



characteristic of a slab standing in free-space. In fact, it 
is impossible to place a PEC plane at the mid-plane of 
a crossed wires slab standing in free-space without dis- 
turbing the fields, because the crossed wire mesh is not 
transformed into itself after reflection with respect to the 
xoy-plane. 



Here, we discuss the application of the proposed ABCs 
to the case where the plane of incidence is the xoz-plane 
(Fig. 2b). The geometry of the double- wire medium 
slab is as explained in Sect. IV. It is assumed that the 
incoming plane wave has parallel polarization (i.e. the 
electric field is in the a;o2:-plane, whereas the incident 
magnetic field is along the y-direction) . 



We have used CST Microwave Studio, in order to par- 
tially validate the dispersion characteristics depicted Fig. 
5, obtained using the homogenization model. Assuming 
plane wave incidence, we have calculated the amplitude 
of the reflection coefficient as a function of the trans- 
verse wave number of the incoming wave ky, for several 
fixed frequencies (Fig. 7). For ky < lu/c the incoming 
wave is a propagating wave {ky = w/csin^^), and, due to 
the conservation of energy, the reflection coefficient am- 
plitude is unity. On the other hand, for ky > w/c the 
incoming wave is an evanescent mode, and thus (since 
an evanescent wave does not carry power), \p\ may be 
any nonnegative value, even larger than unity. Indeed, 
it should be clear from Eq. (31), that the points ky as- 
sociated with the dispersion characteristic of the guided 
modes correspond to the poles of the reflection coefficient. 
As seen in Fig. 7 there is good agreement between the 
response to evanescent waves predicted by the homoge- 
nization model (solid lines), and the actual response cal- 
culated using CST Microwave Studio (discrete symbols). 
In particular, the position of the poles is predicted with 
good accuracy, and thus these results partially validate 
the dispersion characteristic depicted in Fig. 5 associated 
with a = L/10 (dashed line). 



A. Slab standing in free-space 

Considering the geometry of the system under-study 
and the polarization of the incoming wave, it is obvious 
that the magnetic field in all space has a single Carte- 
sian component, Hy. As in Sect. VA, it is assumed 
that inside the metamaterial slab the field is written in 
terms of a superposition of plane waves. The propagation 
constants along z of the plane waves can be calculated 
by solving the dispersion characteristic (26) with respect 
to kz, for a fixed frequency lu, and a fixed transverse 
wave number of the incident wave kx = co / c sin 9i. It can 
be verified that the dispersion characteristic is equiva- 
lent to a polynomial equation of degree three in the vari- 
able fc^ [28]. Thus, the metamaterial slab supports three 
independent plane waves with magnetic field along the 
^/-direction. Clearly, such property is a consequence of 
the nonlocal effects, since in a local material there can- 
not exist two distinct plane waves (i.e. associated with 
different fc^) with the same magnetic field orientation. 
The propagation constants along z of the plane waves in 
the bulk metamaterial are denoted by kz^(i), ^2,(2)1 and 
^z.(3)- Thus, the magnetic field inside the metamaterial 
slab (— L < z < 0) can be written as (the x-variation of 
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the fields is suppressed), 

^-g-'fe.,(2)2 _,_^+g+»fe..(3)Z +^-g-»fc.,(3)^^ 

-L<z<Q (33) 

where are the complex amplitudes of the excited 
waves. For frequencies much lower than the plasma fre- 
quency of the wire medium, uj/c <^ l3p/^/eh, it can be 
verified that in the PEC case and for a propagating in- 
coming wave, only a single propagation constant, let's 
say kz,{i) , is real valued. The corresponding plane wave is 
associated with the hyperbolic isofrequency contours de- 
picted in Fig. 2b. The other two propagation constants, 
kz.(2) and kz.{d,) are pure imaginary, and are associated 
with evanescent modes. 

Assuming that the incoming plane wave propagates in 
the semi-space z > 0, the magnetic field in the air regions 
verifies, 

Hy - H;"=(e^«"-pe-^°"), z>0 

Hy ^ i^;""Te^«^ z<-L (34) 

where ™^ is the complex amplitude of the incident mag- 
netic field, 7o = —iyJuj'^eQ^o — fc^ = yjk"^ — w^eoMo is the 
propagation constant of the incoming wave along the z- 
direction, and p and T are the reflection and transmission 
coefficients, respectively. 

The electric field inside the metamaterial slab can 
be easily obtained from Eq. (33), by noting that 
for each individual plane wave of the type H = 
jij^gife^aigife^z^^^ the corresponding electric field is E = 

HqT^- - I^u^) e»'=-^e*'=-^ where (w, k) is 
the inverse of the dielectric function of the crossed wire 
mesh, defined by Eq. (9) with N = 2. The averaged 
macroscopic current Jd.av inside the metamaterial is fi- 
nally obtained by replacing the electric field in formula 
(6), taking into account that for the present geometry 
k|| = kxiix- The formulas for E and Jrf,av are too long to 
show here, and thus are omitted. 

As in section V, the reflection and transmission coef- 
ficients can be calculated by imposing the continuity of 
the tangential electromagnetic fields {Hy and com- 
ponents of the fields), and the two additional boundary 
conditions (7) {N = 2). Unlike in section V, the two 
ABCs are not degenerate for the present geometry, and 
yield two independent equations. This property is con- 
sistent with the fact that the wire medium supports two 
additional waves for the present configuration, as is man- 
ifest from Eq. (33), and thus two ABCs are required to 
remove the two additional degrees of freedom. 

We have applied the proposed homogenization proce- 
dure to characterize the response of a wire medium slab 
under plane wave incidence. In Fig. 8 we depict the 
amplitude of the transmission coefficient, for a fixed fre- 
quency and lattice constant a, with uja/c = 0.6, and for 
different values of the slab thickness L, as a function 
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FIG. 8: (Color online) Amplitude of the transmission coef- 
ficient as a function of normalized kx for a fixed frequency 
and different L and Eh- The lattice constant is such that 
Lja/c — 0.6, the radius of the wires is r„ = 0.05a. Solid lines: 
Homogenization model. Dashed lines: full wave simulations. 
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FIG. 9: (Color online) Phase of the transmission coefficient as 
a function of normalized kx for a fixed frequency and different 
L and Eh- The lattice constant is such that uja/c = 0.6, the 
radius of the wires is r^, = 0.05a. Solid lines: Homogenization 
model. Dashed lines: full wave simulations. 



of the transverse wave number k^ of the plane wave. 
As mentioned before k^ = uj/csmOi for a propagat- 
ing plane wave, whereas k^ > uj/c when the incoming 
wave is evanescent. For simplicity the wires are as- 
sumed perfectly conducting. The results of Fig. 8 re- 
veal a good agreement between the homogenization re- 
sults (solid lines) and the full wave simulations obtained 
with CST Microwave Studio (dashed lines), both for the 
propagating spectrum (k^ < i-^/c), and for the evanescent 
spectrum (k^ > uj/c). The homogenization results are es- 
pecially accurate in the example where the permittivity 
of the host material is = 2.2. It is seen that the trans- 
mission coefficient may have several poles which indicate 
the presence of guided modes. 

As reported in Ref. [28] , as a consequence of the hyper- 
bolic shaped isofrequency contours associated with the 
propagating mode (Fig. 2b), the Poynting vector suf- 
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fers negative refraction at the interface between a crossed 
wire mesh and air. Therefore, unUke in a conventional 
dielectric slab, the incident beam suffers a negative lat- 
eral shift A as it propagates through the wire medium 
slab [28] (the inset in Fig. 9 represents the case in which 
A is positive). As demonstrated in Ref. [28], the spa- 
tial shift A can be related to the phase 4> = — arg T of 
the transmission coefficient (T = |T| e"**^), by the simple 
formula A = d(j)/dkx. Thus, the spatial shift is propor- 
tional to the slope of the phase of T. Hence, when is an 
increasing function of the lateral shift is positive and 
the beam is positively refracted (this is what happens 
in a conventional dielectric slab). On the other hand, 
when (p is a decreasing function of the lateral shift 
is negative, which indicates that the beam is negatively 
refracted [28]. 

In order to show that our homogenization model pre- 
dicts the emergence of negative refraction, we have cal- 
culated the phase of the transmission coefhcient as a 
function of kx = u;/csin6'i, for the normalized frequency 
Loa/c = 0.6 (w and a are fixed). The obtained results 
are represented in Fig. 9 for different thicknesses of the 
metamaterial slab. The discrete symbols were calculated 
using CST Microwave Studio. It is seen in Fig. 9 that, 
consistent with the results reported in [28], the phase 
(j) is indeed a decreasing function of k^- In particular, 
it is seen that for thicker slabs the slope of becomes 
more negative, which indicates, as could be expected, 
that the lateral spatial shift becomes larger. Similarly, 
the slope also increases with the angle of incidence (i.e. 
with kx), consistent with the fact that for larger values 
of 9i the beam is more refracted at the interface [28] . As 
demonstrated in Ref. [28] , the phenomenon of negative 
refraction is very broadband. 



the behavior of a perfect magnetic conductor. It should 
be noted that in the present configuration the incoming 
wave is TM-polarized, whereas in the Sect. VB the wave 
is TE-polarized. Thus, the high impedance property is 
independent of the polarization. The homogenization re- 
sults concur well with the full wave simulations computed 
with CST Microwave Studio (discrete symbols), which 
further validates our homogenization theory. 
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FIG. 10: (Color online) Phase of the reflection coefRcient as a 
function of normalized the normalized frequency for different 
angles of incidence. The lattice constant is a = L/10 and the 
radius of the wires is = 0.05a, where L is the thickness 
of the grounded slab. Solid lines: Homogenization model. 
Star and triangle shaped symbols: full wave simulations for 
9 — 15° and 6 = 85°, respectively. 



VII. CONCLUSION 



B. Grounded Slab 

For the sake of completeness, we have also studied the 
application of the ABCs to the case where the bottom 
face of the metamaterial slab is grounded {z = —L). As 
in Sect. VIA, the magnetic field inside the wire medium 
is written as in Eq. (33). The magnetic field in the 
air region (z > 0) is defined in the same manner as in 
Eq. (34). For this configuration, the homogenization 
procedure is similar to that described in Sect. VIA, ex- 
cept that at the ground plane it must be enforced that 
Ex=0 and the two ABCs (8) (iV = 2). We have apphed 
such analytical formalism to characterize the response 
of a grounded metamaterial slab with a fixed thickness 
L and lattice constant a = L/10, as a function of fre- 
quency. In Fig. 10 we depict the calculated reflection 
coefficient phase (solid lines) for different angles of inci- 
dence. The host material is air and the wires are per- 
fectly conducting. Similar to the results of Sect. VB, it 
is seen that at certain frequencies, nearly independent of 
the angle of incidence, the metamaterial may behave as 
a high impedance ground plane, and mimic very closely 



We extended our previous work on additional bound- 
ary conditions [35, 36] to the case of wire media formed by 
N different nonconnected components {N = 1,2,3). Us- 
ing simple physical arguments, it was demonstrated that 
in general N different ABCs must be considered at an 
interface. These ABCs were derived for the cases where 
the wire medium is either adjacent to a dielectric or to a 
conducting material. It was proven that in the absence 
of loss the proposed ABCs ensure the conservation of the 
power flow, and in particular we have derived a general 
formula for the Poynting vector in the wire medium. We 
have illustrated the application of the proposed ABCs 
when the metamaterial is formed by a double array of 
metallic wires {N — 2). It was demonstrated that the 
proposed homogenization concepts enable the accurate 
numerical modeling of the considered nonlocal materials 
in scattering problems (for both propagating and evanes- 
cent incident waves) and in propagation problems (cal- 
culation of the guided modes). In particular, we have 
highlighted the anomalous physical properties of noncon- 
nected wire media, and demonstrated that these mate- 
rials may enable the realization of very compact devices 
and the emergence of negative refraction. 
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APPENDIX A: 

In this Appendix wc demonstrate that in nonconnected 
wire media the bulk macroscopic fields may be identified 
to a good approximation with the TA-fields, provided 
none of the imit vectors u„ lies in the transverse (xoy) 
plane. This result is a generalization to double and triple 
wire media of a similar property characteristic of arrays 
of parallel wires [36] . 

To this end, consider the unbounded periodic mate- 
rial and suppose first that the microscopic fields (e, b) 
have the Floquet-Bloch property, i.e. (e, b) e~^^ '^ has 
the same translational symmetry as the lattice, where k 
is the associated wave vector. Following Refs. [18, 37] 
the bulk electric field is defined by E = Eaye*'' '" with 
(using the time convention e~"^*): 

Eav = 77^ / e(r)e-^'''-d3r (Al) 

Vccll J 

n 

The macroscopic induction field Bav is defined similarly. 
In the above, fl represents the unit cell of the noncon- 
nected wire medium, and Vcoii = o,^ is its volume. The 
bulk macroscopic fields verify the system [18, 37]: 

ik X = -iw£oe/iEav - «wPav 
ik X Eav = iiO Bav (A2) 
where the generalized polarization vector, Pav, is given 

by, 

Pav= . t / Jc(r)e-^''-'^ds (A3) 

dD 

and Jc is the density of current over the surface dD of 
the metallic wires enclosed in the unit cell. The dielectric 
function of the bulk material (defined by Eq. (9)) is such 

that f (w, k) .Eav = eoShEav + Pav [18]. 

The objective here is to relate the bulk medium fields 
with the TA- fields defined by (1). Notice that unlike 
the bulk medium fields, the definition of the TA- fields 
(lepends on the considered transverse plane, or equiva- 
lently, depends on the orientation of the wires relative 
to the transverse plane (which is assumed in this work 
to be the xoy plane). We suppose that the transverse 
unit cell Qt niay be related to the unit cell of the bulk 
material as, f2 = Qt x [— flz/2, az/2] for some which 
may depend on the orientation of the wires relative to the 
transverse plane. From Eq. (26) of Ref.[37] it is known 



that the generalized polarization vector may be written 
as a function of the TA-density of current as, 

Pav = — ^ / 3d,.y{z)e-"'^'dz (A4) 

-az/2 

where kz is the 2 component of the wave vector. But, 
as demonstrated in section II, within the thin wire ap- 
proximation, the TA-current is expressed in terms of the 
microscopic cmrcnts as in Eq. (5). Moreover, since we as- 
sume here that the microscopic fields verify the Floquet- 
Bloch condition (along three independent directions of 
space) it is clear that the microscopic currents must be 
such that In{z) = InG^^'^, for some /„. Taking this prop- 
erty into account, and using Eq. (5), it is immediate that 
Eq. (A4) implies that: 

Jd,av(^) = -^^'Pave^'^' (A5) 

Therefore, within the thin wire approximation, and in- 
dependent of the orientation of the wires relative to the 
interface, the TA-current associated with a Floquet mode 
of the unbounded material can be written in terms of the 
polarization vector of the bulk material as in Eq. (A5). 
Such result readily implies that thc^ TA-fields must verify: 

Eav,T(^)=Eave^'^^ Bav.T (2) = Bav 6^ ^ (A6) 

Indeed, in the considered scenario, the TA-fields are uni- 
vocally determined by the (fc^-Floquet) solution of the 
differential system (2) with Jd.av given by Eq. (A5). 
But since the bulk medium fields verify Eqs. (A2), 
and noting that the wave vector can be decomposed as 
k = k|| + kzUz, being k|| the projection of the wave vec- 
tor into the transverse plane, it readily follows that, in- 
deed, the solution of the differential system (2) is given 
by (A6). This demonstrates that the TA-fields can be 
identified with the bulk medium fields, as we wanted to 
show. It should also be clear that Eq. (A6) can be im- 
mediately generalized to the case where the microscopic 
fields are a superposition of several Floquet-Bloch modes 
(associated with different wave vectors). In such case, the 
TA-fields are obviously given by a superposition of plane 
waves (being each plane wave associated with a different 
Floquet mode as in Eq. (A6)). 

As a final remark, we note that the result (A6) is 
only valid provided there are no wires parallel (or quasi- 
parallel) to the interface, because, as noted in section 
II, Eq. (3) becomes singular in such circumstances. In 
fact, if some of wires are parallel to the interface, it was 
demonstrated in Ref.[37] that the TA-fields cannot be 
really identified with the bulk macroscopic fields. The 
reader is referred to Ref.[37] for more details about the 
methods that can be used to study such configurations. 

APPENDIX B: 

In this Appendix we demonstrate that for the super- 
position of plane waves (22), Sz defined by Eq. (19) may 
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be written as in Eq. (23). 

To begin with, we note that the polarization vector 
associated with the superposition of plane waves is: 

P(r) = eoJ2{^i^M - Shfj .E;e+''^'-'- (Bl) 
I 

where e is the dielectric function of the wire medium. In 
particular, the projection of the P onto the direction u„ 
(parallel to one of the wire arrays) verifies, 

Pn=eoY^ i^n,n (w, k,) - Eh) u„.Eie+' (B2) 

and, 

(*ll + ^"^) = ^o'^{i^i-Un) X 

I 

{£n,n (W, k;) - Eh) U„.Eie+''''-'" 

(B3) 

I 



On the other hand, explicit calculations show that: 



{Sn,n - £ft) ^ = (en,n - £ft) (-2k.U„) (u^.Un) 

(B4) 

Substituting Eqs. (B2)-(B4) and Eq. (22) into Eq. (19), 
it is found that: 



Sz — — 



l.m 



Re 



' N 

E 

i^n— 1 



(£„,„ (w.ki) - Sh) 



Bf 1 
(£„,„ (o;, k„) - ShT E^-fln^ (c^, kO u„.Eie^('''-''-)--- \ 



(B5) 



This result and the definition of (11), show that Sz 
may be written as a linear combination of exponentials 
of the type e'^'*'"''"'^ '', where I and m index a generic 
plane wave from the considered set. However, from Eq. 
(18), we know that in the lossless case must be a 
constant, i.e. independent of r. This means that in 
Eq. (B5) the coefficients associated with the exponen- 



tial e''^'^' ^rn)-^ necessarily vanish when k; ^ k^. Hence, 
using the properties (cj, k„j)* = e„_„(a',k^) and 

~ de 

M~ ~ S g£'" UkUk, we finally conclude that may 

n=l 

be written as in Eq. (23). 
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